Nonequilibrium Fluctuation Relation for Sheared Micellar Gel in a Jammed State 
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We show that the shear rate at a fixed shear stress in a micellar gel in a jammed state exhibits 
large fluctuations, showing positive and negative values, with the mean shear rate being positive. 
The resulting probability distribution functions (PDF's) of the global power flux to the system 
vary from Gaussian to non-Gaussian, depending on the driving stress and in all cases show similar 
symmetry properties as predicted by Gallavotti-Cohen steady state fluctuation relation. The fluctu- 
ation relation allows us to determine an effective temperature related to the structural constraints 
of the jammed state. We have measured the stress dependence of the effective temperature. Fur- 
ther, experiments reveal that the effective temperature and the standard deviation of the shear rate 
fluctuations increase with the decrease of the system size. 
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For a system driven arbitrarily far from equilibrium 
by a large perturbation, the traditional linear response 
theory and the fluctuation dissipation theorem do not 
apply. However, remarkably strong Fluctuation Theo- 
rems (FT) have been obtained for a variety of driven sys- 
tems arbitrarily far from equilibrium. Motivated by the 
molecular dynamics simulation results on sheared hard 
disks in two dimensions [T], Evans and Searles derived 
a FT for systems going from an equilibrium to a non- 
equilibrium steady state [2 and Gallavotti and Cohen 
derived FT for non-equilibrium stationary state systems 
[5] . The Gallavotti- Cohen Steady State Fluctuation Re- 
lation (SSFR) based on chaotic hypothesis says [5], 

Ltr^^P{ + Sr)/P{~Sr) = 6^ "^ (1) 
rt+T 

with S = 1. Here, Sr = \ I s{t )dt and s{t) is the 

rate of entropy production in the non-equilibrium steady 
state. P{+Sr) is the probability of observing a fluctua- 
tion of magnitude Sr over a phase space trajectory of du- 
ration T which is larger than any microscopic time scale 
of the system. Naturally, P{—Sr) gives the probability 
of transient violation of second law of thermodynamics 
for the time r, as the entropy decreases over this time. 
The physical implication of Eq(l) is that, if the value of 
Sr and T is large, as in case of macroscopic systems and 
time scales, P{+Sr) » P{—Sr), i.e. the probability of 
observing entropy increasing fluctuations are overwhelm- 
ingly large compared to those in which entropy decreases. 
Thus, in classical thermodynamics we never see the de- 
crease in entropy in any physical process. Extension of 
steady state fluctuation theorem for finite times is dis- 
cussed in [B]. 

The experiments on Fluctuation Relation (FR) re- 
ported so far can be broadly divided into two classes. 
Experiments on systems with small number of degrees of 
freedom include dragging of a Brownian particle in an 
optical trap [3 [8] , electrical circuits [9] , RNA stretching 
[TCTl ITT| where RNA free energy between folded and un- 



folded states were estimated using Crook's Relation and 
Jarzynski Equality and stochastic harmonic oscillators 
[12 [12] ■ The second class of experiments which is of rele- 
vance here, includes macroscopic systems with large num- 
ber of degrees of freedom such as, Rayleigh-Benard con- 
vection [TH [TS] , pressure fluctuations on a surface kept 
in turbulent flows [3], vertically shaken granular beads 
f^, Lagrangian turbulence on a free surface [16 and liq- 
uid crystal electro-convection |17j . To our knowledge no 
experimental evidence exists for the FR in large volume 
sheared fluids. In this Letter we address, for the first 
time, instance of the FR in case of a macroscopic sized 
sheared micellar gel in a jammed state. In our context 
s{t) = -jT^^^i where P(t) is the instantaneous power 
flux into the system and T^ff is the effective tempera- 
ture of the system. We show that the nature of PDF's of 
global power flux for the same system can be Gaussian 
or non-Gaussian, depending on the applied stress [15] . 
Further, the PDF's show similar symmetry properties 
as predicted in Eq.(l), even in non -trivial cases where 
PDF's show large deviation from the Gaussian nature. 
An important point is that in our sheared system, the 
dominant noise is not thermal but rather athermal, i.e it 
arises in the very process of driving, as in [5]. Our results 
also provide a method to estimate the effective tempera- 
ture Tef f of the driven jammed state, like in macroscopic 
granular system [5]. It should be emphasized that T^ff 
is an effective structural temperature estimated from the 
non-equilibrium fluctuations resulting from the applied 
stress and has nothing to do with the actual tempera- 
ture of the system. We will show that Te// increases 
with the increasing driving stress as well as decreasing 
system size. This result may have important signiflcance 
in statistical formulation of driven complex systems. 

Our experiments are carried out on surfactant 
Cetyltrimethylammonium Tosylate(CTAT)-water sys- 
tem in the concentration range 35 - 41 wt% [19] where 
a hexagonal phase of surfactant cylindrical micelles is 
formed. The samples are prepared by dissolving known 
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amount of CTAT (Sigma Aldrich) in double distilled wa- 
ter and are kept for equilibration for a week at 
All the experiments are carried out at 30°C using MCR 
300 stress controlled rheometer (Anton Paar, Germany) 
which has a minimum angular resolution of 0.01 ^rad. A 
humidity chamber is used to minimize the evaporation of 
water from the sample during the rheology experiments. 
We have used cone and plate (CP) and parallel plate 
(PP) geometries with the rheometer. For both of them 
the bottom plate is fixed and top plate (cone with angle 
1.99° in case of CP) rotates to apply shear on the sample 
kept between the plates. For CP the sample volume is 
1.3 X 10~^m'^ and for PP the sample volume was varied 
between 0.98 x lO-'^m^ and 0.49 x 10"^ m'^ as the gap 
is changed from 200 /im to 10 fim, respectively. The re- 
sults presented in Fig.l and Fig. 2 are done in CP and 
the system size dependence (Fig. 3) is done in PP. The 
sample is subjected to a constant shear stress (cr) and 
the shear rate (7) is measured as a function of time. The 
time resolution is 40 ms between two consecutive data 
points in the shear rate measurements. The shear rate 
shows interesting time dependence, which is known as 
'aging' and 'shear rejuvenation' [20]. Here, below a crit- 
ical shear stress, the shear rate decreases with time (i.e. 
viscosity increases with time, known as 'aging') as shown 
in Fig. la, whereas at higher stresses, the viscosity de- 
creases with time (known as 'rejuvenation' which is not 
shown in Fig.l). To characterize the jammed state we 
carried out a stress sweep on the sample (Fig.lc) with a 
waiting time of 20s for each data point, after the steady 
jammed state is fully reached (Fig. lb). We see that below 
a stress of ^--^5 Pa the shear rate is very small but starts 
to increase beyond 5 Pa. Thus, below 5 Pa stress the 
sample is in a jammed state which acts like a soft solid 
under small perturbations. The solid line in Fig.lc marks 
the approximate boundary separating the jammed state 
and rejuvenated state. We also carried out a frequency 
sweep measurement on the jammed steady state as shown 
in Fig. Id by applying a sinusoidal stress of amplitude 1 
Pa. Over the entire frequency range, elastic modulus C 
remains larger than viscous modulus G" signifying the 
solid-like behaviour of the jammed state. The crossover 
frequency is <0.008 rad/s, implying that the relaxation 
time is >120s. The main focus of this paper is in the re- 
gion of jamming for CTAT 39 wt%. For an applied stress 
of 2 Pa, the sample goes into a stress-induced jamming 
phase after initial aging for ^ 60s, (Fig. la). The nature of 
the shear rate fluctuations when the jammed state is fully 
reached is shown in Fig. lb, which shows almost equal 
number of positive and negative values. In the steady 
state the average shear rate is <7>= 3.66 x lO^^s^^. 
The range of fluctuations are much higher than the in- 
strument resolution plj . The global power flux at time t 
from the rheometer drive to the system, P{t) — cr 7 {t)Vs, 
where Vg is the volume of the sample. We define a nor- 
malized variable, W{t) = s{t)/ < s{t) >=7 / <7> and 



rt+T 

Wr = I j W{t)dt . Obviously, Wr = s^l < s{t) > 

where < s{t) >= a <7> Vs/ksTeff and <7> is the 
time averaged value of shear rate in the steady state over 
entire run time. Taking E = 1 as in [5l, Eq(l) can be 
simply written in terms of normalized variable, in the 
large r limit as, 

- ^ — < '«» 

(2) 

Using time-series analysis, we see that the data 
(Fig. lb) do not correspond to low-dimensional chaos [55]. 
From these fiuctuations in the shear rate, we construct 
Wr time series. At this point we want to make a brief 
comment on the averaging procedure. In constructing 
the Wr time series, we divided the s{t) series into dif- 
ferent bins of length r. To improve statistical accuracy 
we have also taken overlapping bins. To ensure inde- 
pendent sampling, the centre of each bin is shifted from 
the previous one by a time difference (0.12s) larger com- 
pared to the correlation time which is < 0.04s. The PDF 
is the same as obtained by non-overlapping bins, but 
with higher statistical accuracy. The PDF's for Wr are 
strongly non-Gaussian for integrated power fiux for all 
r's. We have shown the PDF's for r = 0.4s (Fig.le) 
and T = 0.56s (Fig. If) and the solid lines correspond 
to Gaussian fits to the data. In both the cases, PDF's 
show strong deviations from Gaussian nature for \Wr \ > 
5, but the quantity lii[P{+Wr) / P{—Wr)] goes hnearly 
with Wr, up to Wr ~ 9 as shown in Fig.lg. The result is 
non-trivial and correspond to almost one order of mag- 
nitude variations in PDF's. Fig.lh shows a plot of R vs 
Wr where all the curves scale into a single master curve; 
a straight line passing through the origin, thus agreeing 
with Eq.(2). The corresponding slope gives < s{t) > = 
0.8± 0.006 which corresponds to an effective temper- 
ature T^ff = (8.8 ± 0.07) X 10"K. 

We will now present the result for CTAT 39 wt%, 
when the applied stress is reduced to 0.5 Pa as shown 
in Fig. 2. The experiments are done on fresh samples 
from the same batch of CTAT (39wt%). In this case the 
sample goes to a stress-induced jamming state as soon 
as the experiment is started. The average shear rate 
for the steady jammed state is <7>= 3.77 x lO^^s"^ 
which within experimental errors is same as for 2 Pa 
stress. The typical nature of shear rate fluctuations in 
the jammed state are shown in Fig. 2a. In this case, the 
probability distribution functions (PDF's), P{Wr) of WV 
(Fig.2b) are perfectly Gaussian for all r's (only three are 
shown in the figure for clarity) , as depicted by the Gaus- 
sian fits in the figure almost over three orders of magni- 
tude. From these PDF's we again estimate the quantity 
R for different r's and plot them against Wr, as shown 
in Fig. 2c. All the curves again scale into a straight line 
passing through origin. The magnitude of the slope of 
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FIG. 1: (a) The initial aging part for CTAT 39wt% for the 
applied stress cr = 2 Pa. (b) Typical shear-rate fluctuations 
for CTAT 39wt% for a = 2 Pa. (c) Average shear rate (<7>) 
vs stress (a) and (d) Storage (G') and loss (G") modulus vs 
angular frequency for CTAT 39 wt% sample (stress amplitude 
— 1 Pa ). The dashed lines are guide to the eye. Probability 
distribution functions of Wr for (e) t — 0.4s and (f ) r — 0.56s, 
which deviate strongly from Gaussian nature as seen from the 
Gaussian fits (sohd lines), (g) Plot of ln[P{+Wr)/P{-Wr)] 
vs Wt for different r's: 0.4s, 0.56s, 0.76s against Wt, solid 
lines are the straight line fits to the data, (h) Plot of R = 
^.ln[P{+Wr)/P{-W^)] vs Wr for different r's: 0.4s, 0.56s, 
0.76s. Here, all the curves scale into a straight line passing 
through the origin, as shown by the fitted solid line. 



the straight line gives < s{t) > = 1.2 ± 0.0032 which 
gives Teff = (1.5 ± 0.004) x 10"K. We have determined 
Tef f for eight different values of the applied shear stress 
(cr), shown in log-linear plot, Fig. 2d. It can be seen that 
Teff increases with a, similar to the variation of effec- 
tive temperature with respect to shear rate in We 
have done experiments for a few more concentrations of 
CTAT. For samples with a concentrations of 35 wt%, the 
nature of the PDF's for Wj- remains non-Gaussian for a 
stress level of more than 0.5 Pa. The non-Gaussian fluc- 
tuations of Wr are observed at higher stress values as the 
surfactant concentration is increased. 

We will now turn our attention to the system size de- 
pendence of the nature of the observed fluctuations. For 
this purpose, we have used PP as mentioned earlier, be- 
cause in CP there is a limit to study system size depen- 
dence (when the tip of the cone touches the plate, but 
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FIG. 2: (a) Typical shear-rate fluctuations for CTAT 39wt% 
for a = 0.5 Pa. (b) Probability distribution functions P{Wt) 
vs Wt- The solid lines are the fits to the Gaussian function. 

(c) R = ^.ln[P{+Wr)/Pi-Wr)] vs Wr ioT different r's: 0.52s, 
0.6s, 0.68s, 0.76s, 0.84s. All the curves scale into a straight 
line passing through the origin, shown by the fitted solid line. 

(d) Teff vs applied stress a. 



still appreciable amount of sample remains in between) . 
This problem can be overcome by using PP where the 
gap between the plates can be reduced to a arbitrarily 
small value. We are aware that PP is not ideal for rhe- 
ological measurements since in this geometry, the shear 
rate increases linearly with the radial distance from the 
centre of the plate (which is compensated in CP by ad- 
justing the angle of the cone). We will, therefore, get an 
effective value of the shear rate. Since we are interested 
at this point, not in the absolute values of the shear rate 
but rather in the statistical properties of fluctuations in 
7, we'll sidestep this issue. The experiments were done 
for 39wt% sample for the applied stress of 0.5 Pa, with 
seven gap thicknesses varying between 10 ^m and 200 
/im. The nature of fluctuations in 7 are shown in Fig. 3a 
for three values of gap thickness. It can be seen that 
the amplitude of fluctuations increases as the gap is re- 
duced. Consequently, the width of distributions of Wr 
also increases. For each gap, the probability distribution 
function of Wr (only for 0.72s) is shown in Fig. 3b, for 
clarity. The PDF's are Gaussian for all r's in all the three 
cases. This is shown by the fit of the data to the Gaussian 
function P{Wr) = l/\/2^e- (^--<^-»'/2r' where 
< Wr > and are, respectively, the mean and vari- 
ance of the PDF's of Wr- Here, the Gaussian na- 
ture of fluctuations enables us to estimate the ratio of 
probabilities directly, namely, bi[P{+Wr) / P{—Wr)] ~ 
{2 <Wr > /T'^)Wr. Thus, symmetry hke SSFR re- 
quires the quantity l/Wr[ln[P{+Wr)/ Pi-Wr)]] ^ 
(2 < Wr > /r^) to be a linear function of r. We have 
estimated the quantity (2 < Wr > /tT^) for different r's 
from Gaussian flts and plotted them against r, for three 
different gap thicknesses, as shown in Fig. 3c. In all cases 
we get straight lines parallel to r axis (Fig. 3c) in the large 
r limit, as predicted by SSFR. The effective temperature 
Teff, estimated from Eq.(4) for different gap values (L), 
is shown in Fig.3e. In calculating T^ff we need a value of 



4 




1/L (nrrr^) 



be defined experimentally using FR. It will be interest- 
ing to explore the connection between our results and a 
recent theoretical work on defining the temperature of 
a static granular assembly Our observations opens 
up a possibility of formulation of statistical mechanics of 
a driven jammed state in an equivalent stress ensemble. 
Our preliminary experiments on other systems such as 
clay suspensions showing jammed state suggest that the 
non-equilibrium fiuctuations in viscosity, both positive 
and negative, are generic features of a jammed state. We 
hope that our experimental observations will stimulate 
many new experiments and theories on non-equilibrium 
fiuctuations in macroscopic systems. 

AKS acknowledges Council of Scientific and Industrial 
research (CSIR) of Government of India for support. 



FIG. 3: (a) The time dependence of 7 for different gaps L be- 
tween the parallel plates (L = 10 fim, 70 /im and 200 fim) for 
CTAT 39wt%. (b) Probability distribution functions of W-r 
for r = 0.72s for L = 10 ^m, 70 fim and L = 200 fim. The solid 
lines are fits to the Gaussian function, (c) 2 < Wt > /rF^ vs 
T. The dotted lines indicate the asymptotic convergence of the 
function for large r. (d) Standard deviation F of shear rate 
fluctuations as a function of 1/F. The solid line corresponds 
to F oc 1/L. (e) Effective temperature (T^ff) as a function of 
1/L. The error bars are less than the size of the symbol. The 
dotted line is a guide to eye. 



the sample volume Vs which is estimated from the geom- 
etry of the shear cell. In practice, a small volume of the 
sample can protrude out of the shear cell which will add 
error bars in T^ff, especially for small L. This error is 
not possible to estimate accurately. It can be seen that 
Teff increases as L decreases. The dotted line a guide 
to the eye. Further, the standard deviation (F) of the 
shear rate fiuctuations shown in Fig. 3d, increases as F 
decreases. The solid line corresponds to F oc 1/L. 

To conclude, wc have shown the evidence of Fluctu- 
ation Relation in a sheared micellar gel in the jammed 
state. The system under study is a macroscopically large 
system which is expected to display " entropy consuming" 
fiuctuations if the dynamics is determined by long lived 
temporal and spatial correlations UHl HH]- These long 
lived correlations will be present in the jammed state. A 
physical picture of observing negative shear rate fiuctu- 
ations is that a rearrangement in the jammed state can 
result in feeding the elastic energy from the system to the 
rhcometer drive, resulting in negative 7 fiuctuations. Re- 
markably, we could control the nature of PDF (Gaussian 
or non-Gaussion) for the same system with applied stress 
as a tunable parameter. The symmetry properties of the 
PDF's like SSFR is verified in the non-trivial case of non- 
Gaussian fiuctuations with good statistical accuracy, due 
to huge number of large negative fiuctuations. We show 
that an effective temperature of the jammed state can 
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